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ABSTRACT 


A formalism is developed to describe those 
thermodynamic properties of a metallic system which 
are due to the electron-ion interaction. It is 
shown, in second order perturbation theory, that 
the fee energy due to this interaction may be 
expressed as a simple formula involving the electron- 
ion pseudopotentials and the dynamical structure 
factors of the ionic system. 

The result may be applied to the usual range 
of substances which lend themselves to pseudopotential 
treatment, and possesses considerable generality. It 
may be applied to pure metals and alloys (whether solid 
or liquid)» allowing tor full generality of the,ion 
dynamics, and including those alloys in which mass 
differences and lattice distortion effects are important. 

A number of well known but diverse results are 
obtained as special cases of the general formalism, and 
some numerical calculations are performed as illustra- 
tive examples. 

A particular example is the case of Ag-Au alloys. 
The deviation from linearity (in concentration) of the 
low-temperature specific heat coefficient has been the 
object of some interest in recent years, and a mechani- 


sm, the dynamic "interference" part of the electron-ion 
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interaction, is proposed, which seems to adequately 
account for the observed behaviour, both in general 


magnitude and concentration dependence. 
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INTRODUCTION 


§0.1 General Remarks 


In the past 30 years or so, theoretical under- 
standing of the properties of metals and alloys has 
progressed considerably; and one of the areas in this 
field receiving particular attention has been the 
attempt to enumerate and explain those phenomena or 
properties which arise specifically from the interaction 
between the conduction electrons and the constituent 
ions of such substances, 

The present work will be primarily concerned with 
the development of a formalism which can be used to in- 
vestigate those thermodynamic properties of metallic 
systems which arise from the electron-ion interaction, 

The Ent eth will apply to pure metals and alloys 
and will, in its most general cases, be applicable to 
those substances in which mass, size, and force constant 
differences play a Significant role. 

It has been known for a long time that, in inves- 
tigating the electronic properties of metals, considera- 
ble qualitative success can be obtained by treating a 
metalwas.a Collection of stationary <ions interacting, 
through a weak potential, with a "gas" of conduction 


electrons [1]. 
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The success of this "Nearly Free Electron" model 
would, at first sight, seem to be rather anomalous. 

The electron-ion interaction is essentially a Coulomb 
potential; and for an electron close to the nucleus of 
an ion, the attractive force is very large. 

The first major contribution to our understanding 
of this problem came in 1928 when Bloch [2] showed that 
an electron moving in a perfectly periodic potential 
is not ‘subject to any scattering, apart from’elastic 
Bragg reflections. Thus, in the case of a pure crystal- 
line solid, one may argue that if all the ions were 
stationary, and fixed at sites on a perfect lattice, 
then the valence electrons would travel unhampered 
through the resulting periodic structure. 

It would then be the deviations from periodicity 
due to ionic vibrations, vacancies, impurities etc., 
that cause the electrons to be affected by the ions, 
and to have their dynamical properties altered from 
those of a free "gas". In this way, any scattering of 
@lectrons by the ions would be due to the difference 
between the actual ionic potential and the periodic 
potential of a pure static solid. Under many circums- 
tances one could reasonably expect this net interaction 


to be small. 
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If one turns to amorphous or liquid metals, how- 
ever, the situation should be quite different, if one 
were to use the same line of reasoning. Here one does 
not have the initial periodic structure from which to 
calculate deviations in the ionic potential, and there 
is no reason to suppose that the electron-ion interac- 
tion would not be quite large indeed, 

Experimentally one knows, of course, that while 
therepane wmarked changes,,«the electronic properties of 
a liquid metal are not enormously different from those 
ofthe. conresponding solid,.« The resistivity dof »most 
metals, for example, increases upon melting; but in most 
cases Liguia ~ 20 solid or less, at the melting point 
[1]. This indicates that the scattering of electrons 
does not increase in any drastic fashion in a non- 
periodic structure, and that even in liquid metals the 


electron-ion interaction is rather weak, in most cases, 


80.2 Pseudopotentials 

The resolution of these difficulties arose ulti- 
mately out of the "Orthogonalized Plane Wave" (OPW) 
approach to energy band calculations developed by 
Hérring mSl940' PSy4Pry By 1959)" Privirps ‘and’ Kleinman 
[5] and Antonick [6] had recast this procedure into 


what is now called "Pseudopotential Theory"; and the 
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subsequent activity along these lines has led to a 
greatly improved understanding of what goes on in 
metallic substances. See [7] for a recent review. 

This method starts by dividing the electrons 
in a metal into two types, namely, "core" electrons 
and "valence" or "conduction" electrons. The former 
are bound to, and localized around, a tant dia nuclei, 
while the ‘latter are not bound to “particular ions’, and 
are thus able to move through the body of the metal. 

For the limiting case of no electron-ion inter- 
action, the wave function of a conduction electron would 
be a plane wave. Thus one expresses the wave function of 
a conduction electron as a sum of OPW's, which look like 
plane waves in regions away from the nuclei, but are 
constructed so ‘as: ‘to’ be’ jorthogonal to’ all the’*core’states. 
This orthogonality condition must be Beceem since the 
conduction electrons simply occupy higher states (in 
energy) of the same Hamiltonian governing the behaviour 
of the core electrons. 

Having set up the problem in this way, it is then 
possible to rework the Schrodinger equation for the con- 
duction SUSeELGHS tito a form which is reduced to find- 
ing the rapidly converging plane wave solution to a 
Schrédinger equation with a modified potential. This 


"pseudopotential", then, consists of the original one, 
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plus a non-local part involving the wave functions of 
the core electrons. 

the net -resule rs that the strong part’ of the 
original interaction (which is near the nucleus) is 
removed, and we are left with the effective interaction 
between the ions and the conduction electrons. Now this 
potential will be weak, in the sense that it will not 
have any bound states, but will not always be weak in 
an absolute sense. This need not be an obstacle in 
practical calculations, however, since the NFE model 
will be applicable to a given substance if the net 
scattering of a conduction electron from an ion is small. 
The principal restrictions are that multiple scattering 
effects must be negligible, and that the wave-functions 
of core electrons on adjacent ions do not appreciably 
overlap. | 

One important aspect of these developments is 
that one can visualize, especially in those cases where 
non-local effects are minimal, a metallic substance as 
a collection of "pseudo-atoms" which individually cause 
weak scattering of conduction electrons. This concept 
is not restricted to periodic structures, or those which 
are nearly so, and one is thus able to treat liquids 
and disordered substances on much the same basis as 


requiar solids. in thnewiight of this then, the “nearly 
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free" behaviour of metallic electrons is seen to bea 


reasonable state of affairs. 


80.3 Aim of Present Work 


The present work is intended to apply to those 
metals and alloys where the pseudopotentials are such 
that the scattering of conduction electrons may be 
treated adequately in the Born approximation. 

In the same spirit, a perturbative approach is 
used to calculate the contribution to the free energy 
due to the interaction between the electrons and ions 
which constitute the metallic substance under considera- 
tion. It is well known [8] that perturbation theory 
applied to the free energy has the advantage that the 
contributions in each order do not involve the singu- 
larities which may arise when perturbation theory is 
applied to the Hamiltonian; and as such, is more appro- 
priate for systems having degenerate states. 

it is sfound that, to second order) an the. per=- 
turbation, the free energy of a metallic system, aris- 
ing from the electron-ion pseudopotential interaction, 
may be expressed simply in terms of a scattering func- 
tion [(Q,w), which, for a pure substance, is closely 


related to the Van Hove [9] dynamical structure factor 


S(Q,W). 
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Knowing the free energy $¢, one may then calculate 
any of the other thermodynamic properties which are re- 
lated to it by standard expression, For the internal 


energy and specific heat we have, for example: 


=e 
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80.4 Previous Work 


The specitic-heat-or-metatlize™~systems is’ the“sub- 
ject of much interest, and will be treated in applications 
in a later chapter. It would then seem appropriate, at 
this point, to review some of the previous work done in 
this area, using various techniques, in order to help 
put the results of the present work into proper perspec- 
tive. 

In 1954 Buckingham and Schafroth [10] did a free 
energy perturbation calculation, using a simplified 
electron-phonon interaction to approximate the electron- 
ion interaction in a pure metal. Their object was to 
determine how the electronic specific heat differed from 
the free-electron value because of the interaction. 

They ‘’demonstrated"that' the electronic specific 
heat coefficient was increased appreciably at T = 0°K, 
increased somewhat more as temperature was raised, and 
then decreased, so that at high enough temperatures it 


yielded the free electron value, 
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In the years following this work, a number of 
workers have undertaken investigations with a view to 
obtaining a more complete and accurate understanding 
of this problem. 

These efforts have been directed at various 
aspects of the problem, such as: obtaining results 
that are valid when the interaction is not small, 
largely through the use of Green's function techniques; 
using more realistic electron-ion interactions, usually 
including the effects of Umklapp processes; and consi- 
dering the effect of the direct interaction between the 
electrons. 

The works of Migdal [11], Nakajima and Watabe 
[12], Ashcroft and Wilkins [13], and Rice [14] will 
serve to illustrate these developments; and the situation 
as<of about 1968 is reviewed at some length by Wilkins 
in his lecture volume "Observable Many-Body Effects in 
Metals" [15]. 

At low temperatures, Wilkins shows that the ini- 
tial increase of the electron-phonon enhancement is 
principally *duestona T £n T term, aS was originally 
indicated by Buckingham and Schafroth. 

At high temperatures, his calculations imply 
that the change in the specific heat actually becomes 


negative; and thus the free-electron value is approached 


from below roughly as 1/t*. 
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iInnthefcase @fvymetal alloys, of course, the 
Situation is more complex, and the scattering of 
conduction electrons is altered by a number of impor- 
tant] factoxsnoeThesvaridation of, the,ionic, potential 
from site to site, and the localized constrictions or 
dilations of the lattice, due to differences in atomic 
volume, have a direct effect on the scattering of con- 
duction electrons. In addition to this, the potential, 
volume, and mass differences between the individual 
species alter the dynamic properties of the lattice, 
and this in turn causes further changes in electron 
scattering. 

Two aspects of the alloy specific heat problem 
have received particular attention in recent years, 
namely, the changes due to the potential differences, 
and those due to the altered emission and absorption of 
lattice vibrational quanta. Papers by Stern and by 
Haga [17], as well as the references cited in these 
works will give a good indication of the type of effort 
that has been made to date. The present work will look 
at these effects in a unified way, as a special appli- 


cation of the general results, in Chapter 4. 


80.5 Plan of the Thesis 


The thesis is divided into four principal chapters, 


With a short fifth chapter to summarize the results and 
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conclusions of the work. 

In Chapter 1 we will develop the main result of 
Bane yehesis, and .comment on its applicability. Chapter 
2 will be concerned with the dynamical structure 
factors which, along with the electron-ion pseudo- 
potentials, are of central importance to the formalism 
So Chanrer cl, 

In Chapter 3 we will show how one may obtain 
some well known, but diverse results, as special cases 
of the formalism, and will derive a concise expression 
for the low-temperature electronic specific heat changes 
in an alloy, due to the dynamic part of the electron- 
ion interaction. 

Chapter 4 contains the results of some simple 
numerical calculations. The electron-phonon enhancement 
for ‘Na 2s Bal ree). and the behaviour of the "inter- 
ference" and "impurity" contributions to ‘the specific 
heat is investigated for a series of "theoretical" 
alloys whose pseudopotentials are similar to that of 
Na. 

In addition to this, we calculate the electron 
phonon enhancement for Ag and Au, and propose an explan- 
ation for the concentration dependence of the electronic 
specific heat coefficient in Ag-Au alloys, in terms of 
the dynamic "interference" part of the electron-ion 


interaction. 
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CHAPTER 1 


FORMALISM 


§1.1° ‘Operator Expansion 


In this Section we describe the operator expan- 
Sion that will be needed in the subsequent development, 

Consider an operator A, and an arbitrary function 
of this operator which can be expressed as a power series: 
F(A) = » q_ ahs Then take A to be further given by 
= Ne + ely where i is an operator whose eigen-problem 
is assumed to be already solved: Nolo = hy lors and=¢ 11s 
the usual expansion parameter used here to keep track of 
terms, and will be put equal to unity at the end. We 
will want matrix elements of F(A) between eigenstate of 
he. This expansion was. treated by "Schatroth in°1951: (18), 


fe) 
and the result is as follows: 
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if the |a>'s are also eigenstates of Aye we obtain the 
ordinary Taylor expansion. The diagonal terms are seen 


to; be: 
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Now, writing the free energy and Hamiltonian as: 


= Edy re eb, AP ee-ed H = Hy + eH, 


and making the following choices and definitions: 
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§1.,.261 Hamiltonian 


The Hamiltonian for this problem is given by: 


2 
/ Es ~i 
(he, 2) H =H, (R) +) et ) Hee (Ey 7B) 
au ; 


1 
where ry and p; are the position and momentum coordin- 


oe conduction electron, and the symbol R 


ates of the i 
denotes the set of.all ionic coordinates, including 


momenta when necessary. Hy, describes the ionic system 
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in the absence of interaction with the conduction 
electrons. For the general formalism, it is not 
necessary to specify this operator, and thus the 
final results will be applicable to a wide class of 
ionic systems. 

The pseudopotential interaction for an electron 
at position r Msi det ined by ies .action on an. arbitrary 


wave function n(r): 
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where ~,, is the wave function of the c,"" core Pevet or 

the ion at Ri Note that Be is the actual position of 
ee os is ae } 

the j LON, “And=not. Ghe ay Site on a perfect lattice, 

The effects of electronic screening of the ions, 

as well as any direct interactions between the conduc- 

tion electrons are assumed to have been already incor- 


porated into He by a preliminary calculation, 


ne 
It will be found useful to transform this Hamil- 

tonian, using the second quantization scheme for the 

electrons, .ane’ Co combine) the’ “dragenal’ part of the 


interaction term with the kinetic energy of the unper- 


turbed electron gas. Thus, we may make use of the 
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following prescription for a "one-electron" operator 


& [19]: 
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where | k> is the plane wave state of a free electron 
having momentum fik, and ce and Cy are the Fermi creation 
and annihilation operators. The Hamiltonian (1.2) then 


becomes: 


4 ei 
ads ‘et Cy 1C, A (k,k";R) 


it t Jai 
= Aa CyC,E, + JY Cy C,A(k,k';R) 


The prime on the summation indicates that the k = k'! 
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terms are omitted, and the following definitions are 


used: 
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In the general case VEE and thus A(k,k';R), 
depend on both k and kK! separately, due to the non- 
locality of the pseudopotential. In simpler circums- 
tances, where non-local effects may be safely ignored, 
W, becomes a function of q = k-k'! only, and is) just 


equal to the Fourier transform of the local part of the 


interaction, That is, 
Wastin Ksl).8 > Wie = = 
aun ) sia 


En thisecase (hig A) = ) W, (0), which is just a cons- 
tant, equal to ihe average value of the interaction, 
and simply shifts the zero of energy for the electrons. 
We will use this approximation in the numerical calcu- 
lation! in Chapter 4, but for. the general development, 
the full non-local form will be used. 

As might be supposed, the eigenvalue problem 


for Hy is considered to have been already solved, and 


Hy will be regarded as a perturbation. 


eS Expression forirree, Energy 


Now, the eigenstates {|a>} of the unperturbed 


system can be written as direct products of ionic states 
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sg ee, and the many-body states of the electron gas, 
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and with Hy as given by (1.3), we have immediately 


Enact, 
<a|H, | a> = OURS Oe =eOt Te 


Thus the first non-zero correction to oe will come in 
second order, remembering, of course, that ue refers 
to an unperturbed system, the electron gas part of which 
has the "modified" energy levels given by (1.5). 

In the second order term (1.1), the matrix ele- 
ments of the perturbation may be calculated in the 
following manner: 
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Using this result, we then have 
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where Pn is “the probabilitycof_finding-ethe.tens in state 
[n>, and Pio} in the probability of finding the electron 
gas with the set of occupation numbers fo, }. 
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Here R(t) represents the Heisenberg position operators 


of the unperturbed ionic system, and thus we have: 
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where @ indicates that the principal value of the inte- 
gral is to be taken. This term gives exactly the same 
contribution to the free energy as the first term, and 


thus we finally obtain [20], 
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This formula is the central result of this work, 
and will figure prominently in what. is to follow., It 
is applicable to solid or liquid metals and alloys, 
including those whose constituents have different mass 
and volume; and it allows for arbitrary dynamical pro- 
perties of the ionic system. The essential limitation 


is that the electron-ion interaction must be adequately 
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described by local or non-local pseudopotentials of the 
type defined by equation (1,4), 

Any contributions from elastic Bragg scattering 
should be excluded from (1,12), since these are related 
to the original band structure, and thus belong with 
the modified zero-order results. Other than this, (1,12) 
gives all the second order corrections to the free energy, 
and, by deduction, to such properties as the specific 
heat. 

Now, from the form of (1.12), it is apparent that 
[(k,k' ,w) contains all the information that is needed for 
a particular case, and thus the determination of this 
function will be of primary importance. We shall, there- 
fore, discuss IT in some detail, beginning with the work 
in which a similar type of function was just used in the 


' literature, 
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CHAPTER 2 


DYNAMICAL STRUCTURE FACTORS 


§2,.. Introductory ,Remarks 


In 1954, Van Hove [9] examined the problem of 
the inelastic scattering of thermal neutrons in a pure 
substance. For this case, the interaction is well 
represented by a "local" approximation, and will be 
the same for one ion as for any other; that is, RB 
W(q) for all j. If one calculates the cross-section 
for the inelastic scattering of a neutron from a state 
of wave-vector k to a state k', in the Born approxima- 


tion, then the result may be expressed in the form, 
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where S(q,w) is the Van Hove dynamical structure factor. 
This function is the space-time Fourier transform of 
the time-dependent, quantum-mechanical pair distribution 
function for the ions, and may be written as: 
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where See en indicates an ensemble average over the 
unperturbed ionic system. 

Now, the case of conduction electrons in a metal, 
being scattered by the ions, can be treated in the same 
way in the pseudopotential framework, and one obtains 
ae resulti which is identical with (2.14%), except that 
P(g, 41) fe "Leplaced by (ki k™,@) as given=py CLE?). 

The separation of [ into atomic and structural 
factors, as Boyan by (2.2) in the case of a pure subs- 
tance Frin-the “local "Limit ;"1s"indtrcative *of a “more 
general state of affairs, In the case of an arbitrary 
v-component alloy, with non-local effects, we will show 
presently that r(k,k' ,w) may be written as a sum of 
terms similar to (2.2), involving v different pseudo- 


2 dynamical structure factors (DSF's). 


potentials and v 
In most cases, indeed in all cases, if the pseudopoten- 
tials are taken to be real, the number of DSF's may be 
reduced to v(v+1)/2. These functions will be found to 
depend on k and k', only through the combination gi= 
—— while the separate k and k! dependence which 
signifies non-local effects, will be confined to the 
pseudopotentials i, Ola) 2 

Pmaeamportant feature of being able to separate 
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experimental determination of the DSF's. If one knows 
the interaction potentials for neutron scattering in a 
solid say, then using Van Hove's results (or their 
extension to alloys), one should be able to obtain the 
DSF's from an analysis of thermal neutron scattering 
data. 

Once the DSF's have been determined for a par- 
ticular system, they may be then used to study any 
scattering problem in that system. One need only know 
the appropriate matrix elements of the interaction, and 
thus form P(k,k',w). Having done this, one may proceed 
directly to calculate such things as free energy and 
specitic heat corrections from (1.12), or the electrical 
resistivity of an alloy, which may be expressed as a 


different integral of the same scattering function [21]. 


§2.2 a-a' Structure Factors 


The structure factors of an alloy are closely 
connected to the operators which describe the local 
Sree density ofeach constituent: as) a) ‘function of 
time; and one may express the DSF's directly in terms 
of these operators, in a manner similar to the procedure 
of Van Hove, 

Thus, for a v-component alloy, with each type of 


ion labelled by an index a = 1,2,...,V, we may write 
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the local number density of type a as: 


(2.6) pales ad ? oyegr (Sh en) 
: j ee 


where the superscripted (a) indicates that the sum is 
only carried out over this type of atom. The space 


Fourier transform of this operator is then given by 


ig-r 
(244) N, (a,t) = fax e+ 19 n (et) 


ba iq-R5%? (t) 


e 


and using (1.8), we may thus write, 
(2.5)  A(k,k',R(t)) = J W (k,k")N, (a,t) 
ei 


it being assumed that within an alloy of given composi- 
tion, each ion of a particular type carries the same 
self-consistent potential, denoted here by W, (K,K'). 
Now, using (2.5) and (1.7), and suppressing the 
wave-vector dependence of Wy and Ny for the time being, 


we may write I as, 
at Sear * yt 
w= say e <2W AN (0) 2 EO cs Fe : 


Then, upon adding and subtracting terms involving 
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N1, (O)N, (t), we may rewrite [ as, 


4 i * (a) 
(2.6) r(k,k',w) he Wawa Syq1 (Grow) + S004 (g,w)) 


where we have defined the DSF's 
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(2:58 ) sik) (q,w) - |e e€ <n! (g,0)N , (grt) 


is Nr (G/N, (grt) >, : 


We may note that since n, (r,t) is real valued, then 
equation (2,4) implies that N" (q,t) = N, (-d,t). Thus, 
from, (2,7) .and (2.48), we may.see sehat..the Sua! are.real- 


valued and symmetric (with respect to interchange of a 


(a) 


anda), “whitesthe iso, 
O10. 


are anti-symmetric and pure 

imaginary. Thus if the pseudopotentials are real, or 
ae least satisfy W, = IMle*, where y is the same for 
all a, then the contributions from the terms involving 


a are zero. Moreover, as Bhatia [22] has pointed out, 


because of the ensemble average in the definition (2.8), 


(a) 


aa! a 


the symmetry of most systems would require that S 


( (eyp) (8) 2 4 titsB vy el awe Ww 1 Ri "ty a OR), 
_ so aC Saepoars |, wn 


' 
oo ‘Dae 
e' 120 eft baatieb eved ew Samy 
mi _ eel vid Ay 
; : h ue 
~ Pes inh b= zh ; rt : 
ATA? vy uo. py: awe nS S aa = (Wy Phir gig ) 


Nh 


? at ga 1B) id " 


(tap) pi (Oph as? Ot" s ae é op oh 


- ip ca a 
n° DE GREG Dh om 


Ay 
nett ,beulsy [sex Si (ty 2) 8 oer sat acon sae ale 
aaciel , (B= ) Pile (ap) tedt asian Sse oksnupe 
+iései 545 ioe 2d jade Pore ‘yom: ow , By 8) bane: 4, 9) aoe 
» to spashoretal. ot Sgaqace oid bw) obterane bap’ corer | 
ating: Lied diadsimge—E ine Eye sha ad otidw , (Ne 4 

10. ~tgax. exe aistiaszoupboeeq ony At, aur eras paid, 
to} Sits! age ery Stony”, a) UNL = WouReiten sapat a8 | o 
pavtoyal anaes ene. now? aholind byob ifs eds -? the Mig p 
pititoa aad ($5) BidBAe ex .ze¥veet6M 19S, 238 we 
(86S) ‘nodding |b. ade ni “spsveve ecdaGane: ‘Bas 20. panetiod 
0 = (le tead ‘exlipss binow amednye som, Xo eachacl a a 


wo Be: 


30 


At any rate, we shall not include the anti-symmetric 


terms, and thus obtain, 


Gog Wiki wls  ) WO ke 


This method of writing [T is simple and direct, 
Since it gives the de-composition in terms of the corre- 
lations between the positions of the different types of 
ions which make up the alloy. There are other ways of 
constructing IT however, and one may define different, 
but equivalent sets of DSF's, which might have the 
effect of accentuating some physical features of par- 
ticular interest, or which may be more easily related 
to some experimental data. By choosing the proper set, 
it may also be possible to reduce the number of struc- 
ture factors needed, if some parameter, such as the ionic 
volume, for example, is roughly the same for all the 


constituents. 


§2.3 N-C Structure Factors (Binary Alloys) 

For a binary alloy, Bhatia and Thornton [21] 
introduced a set of three DSF's which express the 
dynamics of the alloy in an instructive manner. They 


wrote I. as, 
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(2.10)  Pig,w) = (W)* Soo (gy) + (W,-W,) 75.0 (qu) 
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+ 2 (W,-W,) 5.0 (q,w) 


where W = cw, ae CoWo, and c, is the concentration of 


type a. Sun is quite similar to the single structure 
factor of a pure substance, and describes the correla- 
tions between the fluctuations in particle density. 
Sac describes the correlations between concentration 


fluctuations, and Suc is the cross-correlation between 
the two. 

These DSF's will be properly defined later on, 
but for the moment we may make a few over-simplified 
remarks about the physical significance of the different 
parts of (2.10). The first term represents a "dynamic 
vVirntual crystal” *comtribution, in which the ionic vibra- 
tions are taken into account, while the second term 
arises basically from an "electron-impurity" interac- 
tion. This latter term may be often well represented 
by the “static approximation”, in which its frequency 
dependence is given by a factor 6(w). The third term 
reflects predominantly those effects which owe their 


origins to differences in size and shape between the 


two types of ions. It will be generally negligible if 
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they have close to the same volume, thus leaving 
only two independent DSF's in such cases. 

In a paper subsequent to [21], Bhatia and 
Thornton [23] derived exact expressions for the high- 
temperature, long-wavelength (q v 0), Limit.of the 
"static structure factors (Sry 6G) - [au Say (G7) SEC.) 
from. fluctuation theory... For theucase of a cubic 


crystal, the explicit expressions are: 
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Suc (y) = - SBnk (y) Sag Cy) 


where y is a Un&tevector in the direction of q (remem- 
bering that ea A) se & is the Gibbs free energy, Be 
is the isothermal bulk modulus, Cy is the concentration 


of species 1, and o is the stress tensor. <K«(y) is given 


by the expression 


5 
K(y) = y Lenny) Son. 
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nw 
where se E is the (unit) polarization vector for a 
v 


wave propagating along y with velocity v peand 


yré 
p is the density of the substance. The symbol 6 is 


defined as 


where Toes (€55/9C ) 9 oon , With e¢ being the strain 
tensorsn Lf ene strains: rntroduced bythe -difrering 
volumes of the alloyed species are quite small, then 
6 is closely equal to (Vi -V5)/ (CV +C5V5), where v; is 
the molar volume of species i. Thus one readily sees 


that the requirement for Sy to be zero in this case is 


c 
essentially that Vi= Vor as we have mentioned above, 


We also note that for this case 


CE 24 


and thus for a random alloy S = (Cl lL-C)s. 


eC 
For shorter wavelengths (larger aq), exact expres- 
sions are not available (in explicit form) for these 
structure factors, but as noted in [23], one may obtain 
some qualitative information from the results of X-ray 


scattering theory, which essentially deals with the 


Static structure factors. In particular, one may deduce 
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from the works of Huang [24] and Krivoglaz [25] that 


fon an elastically ssotropic. random alloy, 


where G is the reciprocal lattice vector closest to qe 
and C1 is one of the elastic constants. We may note 
that this expression is again equal to zero for ee yc 
In addition, we see that this result becomes 
infinite fort, = G, and one would thus expect large 
contributions to any quantity involving S for values 


NC 
of g in this region. We will encounter similar behaviour 
in Chapter 4 when we consider the umklapp contributions 
to the electronic specific heat. The principal diffe- 


rence is that for the present case we have (for gq parallel 


tor G)': 


Siele ice ai, 


while in Chapter 4 we will have (loosely speaking) : 
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§2.4 N-C Structure Factors (y-Component Alloys) 


Now, proceeding to the case of a V-component alloy, 
we may follow Bhatia and Thornton, and relate the DSF's 
to the fluctuations in the local number density of the 
alloyed species, rather than the full number density. 


In fact, if one describes the local fluctuations of type 


a by, 
Cneete, cy, am! (rr, ce) = a 
Ga « sees oi 
ider 
N (art) = far 2 n, (r,t) 
(a) 
- ya) Rigi e herbs. onl ts 

ch ae wie) 
‘ a 


where ny = N/V and further recognizes that elastic 
forward scattering (k = k") has no effect on-the 
structure-dependent electron-ion processes, then the 
preceding development, from (2:3) to (2.9) ,, goes. through 
exactly as before, qexcept that there are a lot of terms 


involving 8y , which are simply ignored, 
; v 


ue ok k! 
Now, tm ‘order to -Gonstruct agset fof DSF's which 


S and S in the binary case, we will 
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define an operator which describes the fluctuations in 


reauce”™to s 


local coéncentration of a given speciés by, 


6c, (r,t) = (V/N) (én, (x,t) —c,6n(r,t)) 
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where §6n(r,t) = ) én, (x,t) is the operator describing 
QO 
the total fluctuation in number density, and Ca = N/N. 


The Fourier transform of 6c, (x,t) is given by, 


Lone ‘ 


1 re 
2.11) ce, (qit) = G) fare *™ gc, (ny) 


II 


1 
(y) CN) (qt) es c N(q,t)) 
wien N(q,t) =) N, (art). Also note that, 

Zz z = 
(2 U2) ) c, (r,t) = j) Cc. (q,t) =p Ou, 

a i al 5 
Now, using (2.11), we may write, 
N, (art) = NC (q,t) + c N(q,t) 


and thus, from (2.5), we have, 


A(k,K';R(t))o= W(k;k")N(q,t) +N J W, (krK')C. (qt) 
sl ope er leas z alee 2 


with W(k,k') = ) cW, (k/k"). 
in case of a binary alloy (go. =.1,2),,;we see from 
oes that C, (q,t) = ~C, (q,t) = C(q,t), and that we may 


thus write, 


A(k,k';R(t)) = WN(q,t) + N(W,- W,)C(q,t) . 


In the general case also, it will often be advantageous 
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to work in terms of differences, rather than full 
pseudopotentials, and thus, taking advantage of (2.12), 


we may write 
A(k,k';R(t)) = WN(g,t) + N }(W-W)C, (q,t) 
a. 


where v may represent any one of the alloyed species. 


Using this relation in (1.7), we then obtain, 
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Then, following the same type of argument that preceded 


(2.9), we-may “symmetrize”” this form of [f, and obtain, 
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By expressing Cy (a,t) and N(q,t) in terms of the 

Ny fart), one may easily verify that going from (2.13) 

to (2.14) is exactly equivalent to the progression from 

(2.6) to (2.9), and that one result implies the other. 
Just as there are y(vtl)/2 independent DSF's of 

the Sua! type, there are the same number in the set of 


number-concentration DSF's, of which 


1 is of the (NN) type, 


(v-1) are of the (NC) type, 


v(v-1) 


5 are of the (CC) type. 


One may also, of course, write a given set of 
structure factors as a linear combination of any other 


properly defined set; and for the two types just given, 


we have, 
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and the inverse relations, 
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from which it is apparent that 
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as is required by (2.12) and the definitions (2.15). 
This is simply a statement of the fact that, while 
there are localized fluctuations throughout the subs- 


tance, the total number of ions is fixed. 


62.5 Structure Factors for a ‘Substitutional Alloy 


Although one might hope for progress in the 
experimental determination of the DSF's, it is at 
present a difficult undertaking to attempt accurate 
theoretical calculations of these functions.  Never- 
theless, WeWis Stillspossiblé tolobtain? useful dnsight 
into many physically interesting processes by consider- 


ing somewhat idealized situations. 
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One such example, which we shall make use of 
later on, is that of a substitutional alloy. The 
basic assumption operating in this case is that one 
may interchange ions of the various species at random, 
without causing significant changes in the local ionic 
structure. This implies a certain degree of uniformity 
in size and shape, of course, but is not entirely un- 
realistic. One knows experimentally that it is diffi- 
cult to form homogeneous alloys from constituents which 
differ markedly from one another in this respect, the 
Hume-Rothery 15% rule [26] being the classic indicator. 
This rule states that if the constituents differ from 
one another by more than about 15% in linear dimensions, 
then their mutual solubility will be severely restricted 
(i.e. to less than about DSi. 

in order to introduce "the”“degree of randomness 
required in this case, we take the probability that a 
given position in the alloy is occupied by an ion of 
type a as being equal to ore the concentration of that 


‘species in the alloy. We thus obtain the configurational 


average: 
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This result may be cast into a variety of equivalent 
forms, and for présent purposes, in order to facilitate 
comparison with (2.14), we introduce the pseudopotential 


differences UW aW di and obtain 


* * 
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Say (92%) ¢ as defined here, is exactly the same as (2.15a), 


of course, and by comparing the forms of (2.14) and 


(2.18), we may deduce that 
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which means that the fluctuations in number density 
and concentration are independent of one another, as 
might be expected from the discussion of §2.3 and the 
assumption of no lattice distortion in the present 
case. In addition, we see that 

Saxe (giw) = ey (6yry — Cyr) Sz (Gru) 

ao 
which shows that the concentration fluctuations are 
all correlated through the same dynamical behaviour, 
with each contribution to the scattering being weighted 
by the concentrations of the various species, and by 
the relative strengths of the pseudopotentials. Thus 
in this simplified case one has only two independent 
structure factors to deal with. 

The first and second terms in (2.18) correspond 
exactly to the "coherent" and “incoherent® scattering 
cross-sections that one encounters in the current treat- 
ment of thermal neutron scattering [25]. Indeed, one 
May see trom (2.16) that the coefficient of S,(q,w) is 
just (we - Ww), which is the form more commonly seen in 
neutron scattering theory. The W's are different, of 
course, but the form°-is the same, and the structure 


factors are exactly the same, as has been stated earlier. 
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CHAPTER 3 


APPLICATIONS 


Sao Preliminary Remarks 


For a given system, represented by its appro- 
priate pseudopotentials and DSF's in the scattering 
Sunctiian P(k,k',w), expression (1.12) gives a straight- 
forward prescription for calculating the second order 
corrections to the free energy. A significant aspect 
of this result is that it embodies, in a simple concise 
expression, many results which are usually treated quite 
separately from one another, using different techniques. 
These techniques are often difficult to extend without 
a complete reformulation, while (1.12), within the basic 
limitations discussed earlier, is very general. 

As illustrative examples, we shall consider a 
number of representative cases, which will serve to 


demonstrate the utility of this formula. 


§3.2 Harmonic Solid 


As a first example, let us approximate a pure 
metal by a system of ions executing harmonic vibrations 
about equilibrium positions on a perfect lattice, which 


is the usual "phonon" model. 
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In this case, the Hamiltonian is given by 


SS) 
— 


i . 
Hy = L hw, (asa, + 


where s indicates one of the 3N vibrational modes of 
the crystal. This index may be further written as 

Ss = (Q,&), with each of the N wave-vectors (Q) corres- 
ponding to three polarization modes (& = eg having 
unit vectors a. and frequencies Woe The Bose creation 
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and annihilation operators, as and ass satisfy the 


commutation relation 


1 = 
[a.,acel “ oe 


Then if a is the ee ionic site on the perfect 


lattice, we may write 


Rete.) = Ss. oe Cte 
R(t) = x, + u(t) 


where a is given by the standard theory [27] as 
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with M as the ionic mass. From (3.1) we may readily 


deduce that 
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which is an ordinary c-number. Thus we may write 
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where we have used the common identity ee 
Whachits velisd Ll -[(A,(A;B)]) = [B,[A,;B]) ="0:; 
The exponent of the term to be averaged can be 
thus seen as a simple sum over phonon states, with each 
term commuting with every other term. Then, since we 
may write an ionic state, in this approximation, as 
Tie % 

= 2 apie A 
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average becomes a product over all modes of terms of 
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Using this result in (3.2), we obtain 
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We must then insert this into formula (2.17) in 
order to obtain r(k,k',w). Because of the difficulty 
in performing the time integrations, however, one may 
resort to the procedure of expanding the second factor 
in a power series and performing the integrations term 
by term. Then, making use of the fact that 

te gate a =n? 5 5 

Seals G 
where G is a vector of the reciprocal lattice, we keep 
only first order terms in the expansion and thus obtain 


the “one phonon approximation": 


n © 
a eet ee se ee) on) Occ 
2 Aue 
h(q-é ) “s fi(q-e.) ae 


+ n.6 (whos) 6449,g)} 
Zz ) Woe ; = = 
P [Wik yk") | G,s 2Mw OE ee ae eee Hehe ARSE site 


as 


a Cem, ,x) er, att) t= — j 
| (it mS) : oof + 


; : 
ak (FES) elorrvot atat aiid axgens neds veym av 


YtiW>L2T ib. ett To seussoF olin Mh 1) sausage Od. sas 
VSM sui \cravawod dnebebaeetme omit ont eriieronsg ak 
1osI8% hnoves ods palbnsata 2c Sthbsoory ont bd stoner 
meso Rrioitespetat oddevuaimsoiseg Dees estos yowog Ran . 


ae 
ssid 338% sfit to Sen pokAee emedT 


} 


‘, 

“ - 6 x sy = tr | i 2 {. 

Se i ‘fx 

~~ - = ot “a 
: os 


qeouk gw ,ootsisl Levoserost -on3)> To aie hol nt * ost 
Weide eurtt Sos noLensqxs aS rel emred pone esta ‘eg! o 


er Saeco eno”, a 


i ee 
“ j ; af 
= go> pat yo 7 
240+ piles Ww" “aXe, fe) 4 aS & PANE 


where we have dropped the elastic Bragg terms Soe G 
a 
in line with the remarks made earlier, and used 
a 2 A 2 
~e = e ‘ 
(g-8y ¢) (g-8_g ¢) 


The higher order, or "multi-phonon", terms in 
the expansion of (3.3) may be obtained directly, but 
we will confine ourselves here to first order. One 
point to be noticed, however, is that the "Debye- 
Waller factor" 

n(q-é_)* 


2NMw , (2n.+1) 


-2D (q) 
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= 


need not be expanded for the purpose of performing the 
time integrations, and may thus be left in exponential 
form, as is the practice in X-ray scattering problems. 
But there have been arguments to the effect that this 
term and ae time-dependent terms may substantially 
cancel one another in higher orders [29]. This matter 
remains somewhat unsettled, however, and for present 
purposes we will adopt the policy of expanding (3.3) 
‘Consistently to the same, order, ingthesexponent, treat- 
ing the Debye-Waller part in the same way as the rest. 
Note that using the definition of yey and the 
relation f(a)é(x-a) = £(x)6(x-a), we may rewrite (3.4) 


in the alternate forms 
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which are sometimes useful for calculations. 
Inserting (3.4) into (1.12) and performing the 


frequency integrations gives for the free energy correc- 


tion 
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This is a more complete version of the eee on used 
by Buckingham and Schafroth [10] in their specific heat 
calculations. In their case the Umklapp (GA 0) terms 
were dropped, an isotropic Debye model was used for the 
phonons, and the factor |wik,k") |*(q.8,) °/2Mw was 


replaced by a simplified "interaction parameter" which 


varied linearly with |Q|. 
Expression (3.5) may be separated into two terms, 
one of which depends only on "virtual" phonon processes, 


while the other depends directly on the phonon occupation 


numbers. This gives 
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where it is understood that the sums will eventually 
be converted to integrals and principal values taken 
when the singularities present problems. 

The second term is zero at T = 0, and does not 
contribute to the electronic specific heat coefficient. 
It may be interpreted as representing a renormalization 
of the speed of sound: “The first term in’ (3.6) ‘will 
yield the "electron-phonon enhancement" of the electronic 
specific heat, which is finite at T = 0, and becomes 
zero at high temperatures, as discussed in 80.4. (See 


(10) “and? ibs). 


b2m3 -ebiquids Metal 

Beyond the melting point of a metal, there wiht 
be some degree of short-range order, but there is no 
regular lattice, and such concepts as band structure 
and Bragg reflections are no longer applicable. In 
such circumstances, the scattering function must include 
all elastic’ and inelastic processes. Indeed, the elastic 


processes will tend to dominate the scattering to a 
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large extent, since they may no longer be assigned to 
the band structure. This fact makes it easy to obtain 
a good first approximation to the electronic properties 
of liquid metals. 

In setting up an expression for the electrical 
resistivity Of a metallic subStance, using the current 
variational approach, one finds that the scattering 


function enters in an integral of the form 


foe) 


(erm) I(q) = | au (Bae) T (q,w) 
e - 1 q 


— CO 


in the local pseudopotential approximation [21]. Thus 
at liquid metal temperatures, and for the range of ener- 


gies accessible to the system, 


Bia <<—r14) meng (Gq) l= | dw T(q,w) . 


In the case at hand, we are dealing with the 
same scattering processes that are responsible for the 


resistivity, and in the same approximation, we may use 


S(q,w) = S(q) 6 (w) 


wherers(q) “isthe “static™ structure’ ‘factor, S (q) = 
Jaw S(q,w), which gives the total scattering in a given 


direction q.",Using this relation in (1.12), we obtain 
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and ney) = (v/ (2m) 3) 4nk*/ (aE, /dk) is the density 

of states, for each spin, of he "modified" free 

electron gas, with Pk given by (1.5) as ce 2K */ 210) +P 
Assuming that the liquid metal is still much 

below the "Fermi temperature", we may expand (3.9) in 


the usual manner [1], and obtain for the first two 


terms: 
2 
03. LO} o,=2 dE ny (Ey) G (Ey) pM (k,T) ree UE eB, 
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Thersecond, term anic3.10Jegives a correction 
EOUENe Specific heat, according to the formula AC, = 


2 Z ea : 
-T (0 o5/9T dy , which is 


e€ = t 
Coes) AC / Cy, = -G (uo) 


; e 2a 242 ssid 
with Cy nop ken (uy) tT as the specific heat of the 
free electron gas. In view of the high temperature, 


however, this correction will be completely dominated 
by other effects, and may be safely ignored. 

The first term, however, is of more interest. 
It is seen to be independent of temperature, and since 
the total energy is equal to the free energy at T = 0, 
we may, in view of the form of (3.8), interpret this 
term as a correction to the electron energy. Thus we 
write the energy of an electron in a state k, to second 


order, as 


oo) Jw Kk") |*8(q) 
e(k) = FE, + = = 
is kséeN Kk! (E) Eye) 


~~ 


This is) the usual result that is normally obtained from 
perturbation theory applied to the Hamiltonian. It may 
be used, for example, to calculate the electronic 


density of states in a liquid metal [30]. 
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§3...4 Binary Alloy =- Static Effects 


From expression (2.18), we may write the scat- 
tering function for a binary substitutional alloy, 


with c as the concentration of species 2, as 


(3212), Tikka) = |e.) |? Sq, u) 
+ c(1-c) |W, (k,k") - Wy (k,k")|* S,(q,w) . 


The first term as jhe’ “virtual crystal” contribution 
to the scattering, and its behaviour will be much the 
same as that of a pure substance with an average 
potential W associated with each site. 

The second term represents the scattering due 
to deviations from the mean potential. It is a sum of 
self-correlation terms, one for each ion, and thus 
displays no "coherence" effects. The static part of 
S, (a,) will often dominate the contribution from this 
term, as may be seen from the following considerations. 
If we consider a binary alloy in which the two types 
Of, ions differ’ only in -itheir respective) potentials, 
then in the one-phonon approximation, we obtain from 
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Using this, we obtain, after some rearrangement 


BPs = ae. 

S>(G,w) dw ar h(q-e.) w, (ng+1) won, 
Poe a Te 8 ome hee, Ooo, | 
ik ray S ad Sy Kk uc ie! 
The term involving the phonons may be estimated by 


using an isotropic Debye model, and at low tempera- 


tures, its maximum value will be of the order of 
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In those cases where this estimate is reasonable, we 
may replace S7 (dr) Mise Oy Oi ane: cous Obtain 


for the free energy correction due to this term 
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Then, using the same procedure as was used to obtain 


(3.11), we £2nd the correction to. the, electron: specific 


heat arising. £rom this: term: to be 
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This result was obtained by Stern in 1966 [16], for 
the case of a binary substitutional alloy, in which 
the ions were fixed at sites ona perfect lattice, 

Stern applied this result to the case of Ag-Au 
alloys in an attempt to account for the deviation. of 
the*electronic’ specific heat’ from’ Llinéar interpolation, 
in concentration, between the values for pure Ag and 
pure Au. Subsequent measureme nts of the low tempera- 
ture specific heat by Davis and Rayne [31] indicated 
that while the concentration dependence of this result 
is essentially correct, its size (as estimated by Stern) 
istoosmall, by about a’ ftactorof three? 

The reason for this lack of agreement is to be 
found in the neglect of the ion dynamics. While it 
is true that the term involving S7(q,w) ire (SSL2) “may 
be adequately treated in the static approximation, as 
we have roughly indicated, this does not mean that one 
may treat all aspects of the electron-ion interaction 
in this way. At low temperatures, ‘in fact, the dynamic 
effects, usually referred to in terms *of the electron- 
phonon interaction, are commonly known to be quite 
important, and we shall return to this point after 
obtaining some necessary formulae in the following 
sections. In particular, we shall perform some numeri- 
cal calculations for the Ag-Au system, in order to 
indicate the importance of the electron-phonon effects 


in this case. 
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§3.5 Low Temperature Specific Heat 


In this section we shall derive an approximate 
expression for the correction to the low temperature 
electronic; specific.»heat coefficient “y" due to the 
term involving Sx (Gr) « The unperturbed value a is 
taken to be the free electron value in the case of a 
disordered system, or the zero-order band structure 
value in those cases where this concept has meaning. 
In either case, the Sun term represents scattering of 
electrons around the Fermi surface due to the emission 
and absorption of lattice vibrational quanta, and the 
fact that the Fermi energy is much larger than the 
energy of the vibrational quanta is an essential point 
in the derivation. 

In addition to this, we will use the fact that 
the low temperature specific heat is proportional to 


the density of states, and use the approximation 


An(EL) 
elo pe AE (Kk) ] 
(EX) (0) (Ox 


where A indicates the change in each quantity from the 
unperturbed value. AE (k) will be obtained by noting 
that at T = 0 the free energy is equal to the total 
energy, and using Landau's "quasi-particle" prescrip- 


tion for the energy in a state k: 


56 


de 


wt eas 
a 


nsioasidiaais’ ae avizeh tisietanas ! 
ouusesoqniad Wet eit of cnn a 
oid 6d oub “\" Juabottasos Beed. pina 
at s#Hisy. hodupsssaqan oat “taal eid 
& XO #2869 Sly ms sulev nosztoels | ott, od 
suusoussa baad sobio“os@s sd5 10 cmsage bs 
-piiiisem apd. tqeoioo 23 Aste | Stonw aoanp g2ond2 wk 
to poixetssoc atqsestgsx M99 ye Said .38K9 — 
woibeing sdt ot sub svstive. ingot SAP BaNdus 2 on” 
| of) bre (sine teneeieedey Sobstek $5. rey ll 
sft osdt sepzet doum, ei Yowsns Lite. oth, ele 
PtLOC iniGucive tig 2 Stnaup Larcoftssxdiry aig. oe. Pane 

- stobsey tit eer: 

isait tosi efd sey {liv ew ,erng- OF nolgibbe AL 
oJ fsnoidzoqoxq et tned er ou satsamed ‘ok. 
aoitemixorggs ‘aad seu bas ,esiete Hie yoteneb # 


a SS af ) 
(aaa ° eh 
7 CAAA - oe eS eres? ee 
da . n~ = 


ol 


sft mos? yi tdnsip ieee ni Sene((>: ols vossoknni A 
a yd hentetdo saree (aaa “oulay Bodsuda 
WY 

. Indes: ona ot, isips ai WeTSNs ssit orem = 7 _ 


~gtroaeie- hatoiaxaq- iesmp” = ebrdut ia ne .YpTs 


we 


fila | 
sal inal J 


Part of the development will follow a procedure used 
successfully by Ashcroft and Wilkins [13,15] in an 
examination of the electron-phonon enhancement of the 
specific heat in pure simple metals. 


For the free energy correction due to Suny (G1) 4 


we have 
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and we thus obtain directly, omitting the superscript 


on p(°) for convenience, and taking account of the two 


spin states, 
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Changing the sum over k* to an integral, we have 


where ¢ (EL an E_). Using the 6-functions to define 


k' 


the limits of integration explicitly, and treating 


the first and second terms separately we obtain 
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where it has been assumed that Ee is much greater than 


those values of Ex for which the integrand is appre- 


~ 


ciably different from zero, and that the zero-order 


density of states does not vary greatly over this inter- 


val. 
From this expression we then obtain, 
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where we have neglected any explicit energy dependence 
of the pseudopotentials. Consistent with this, we 
will treat the pseudopotentials in the local or "on- 
Fermi sphere" approximation in the rest of this devel- 
opment, so that they may be written as a function of 


q = k-KkS.omlly. Then, introducing the LUNCtLOn, 


~ 


~todnd eine tSvo ei igete yasv Yon. | dob ebbade 


neds xossore eet at 4a ‘gett mach ; ‘5 
-omqqs oi bnexpesni end aot a we 


ROR std feds bets .O4a8: me eM | 
aa 


maT 


 Tistde, neds aw dotardxque elds mont 


Saar : trate 
Ox, ‘ Witt ~~ 
a) 
¢~ 6 | 
a arena (aD) > 
ita) a ae bh” L 
2" te ) 


. BOttemnsgead yousAas tidilaxs WTS bsvbetpen ambit s ow 
Sw , aids slo tw, dase eae Sibtdosyoqobuean “oe 

| a6" NO"LESOL. oda! 'nt Siessdesoquoupng’ ead gaen3)2 

~loveb nity to Jas Sita nk nots smiomias Noccbtige 


wot sy 
ae nose & Be ical ad vent “yo: 38, . 06 an 
eee ea “pabsuboa tgs othe pen 
ne | i oe 
| a 


2n , (Ep) Sy (GW) 
N 


nig) ee [icq |?) au ae 


we are led to the following results: 
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dQ AQ 
— | k [= F (k-k') 


4q 4t 
dQ, da 
1 2 2 k (OK! 
eT [ax k § (K-Kg) | de"! 6 (k'-kg) |e) aE F (k-k') 
F 
. hang | ak] a°k'F (kK) 8 (Keg) 6 (eK) 
(4nk5) Dap | 
es +a a7 [ak F (q) 8 (k-k,) 8 (|k-q|-k,) - 
(4nk=) 3 thi 


q<2k, 


The integral of the delta-functions over ak may be 
evaluated directly as 


chet co 

4, 
| a6 du [ax k76 (k-k,,) 6 ( (k*+q7-2kqu)? - ke) 
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Using this result in (3.14), separating the angular 


variables, and defining a function 


aQ 
aS 
(3-15) Fay (a) la F (q) 
2n, (ED) an a 2 Sy (0) 
; ed |W(q) | fau saaiae ei 
we obtain the final result 
2k 
Ay, NN oh 3 
(3 16) (—-) = | chi Fe. ta )* 0. 
Y 2 AV 
re) 2kn 0 


The form of this result is seen to be quite 
similar to the variational expression for the resis- 
tivity of the type of substances which we are consider- 
ing here [21]. The differences are that in the 
resistivity formula there is an extra factor q? in 
the integrand, and F (q) is replaced by a function 
proportional to I (q) gavem by *(3.7). 

Indeed, in the present case, one may replace 
lWi(q) |*S. (gro) by the full scattering function I (q,w) 
as long as any "static" contributions, involving $(w), 
are treated separately, such as in the preceding 


section. 
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CHAPTER 4 


NUMERICAL CALCULATIONS 


4.1 Electron-Phonon Enhancement for Na 


In the one-phonon approximation, we may use (3.4) 


to obtain for a pure substance at T = 0: 
Anne 2 
(4.1) fau ——— = —j— 6__ 
ee Ss 2Mw~ q Q,G 
cS 
A 2 
oe aac, a 
0 G,é 2Mw~ 
=" q-Gre 


On the other hand, if we calculate the structure factor 
Siq) = I(q)/|w(q) |? appearing in the resistivity formu- 


la, we have 


stg) = [au p88 stg.0) 
= a hii et ~ 
h(q-6_)° °q-0,6 
= 26h eGR CONT MePPN ETE Foca CREAN OCT, ES 
cheove 28. fox ghfusmangya leet nesy 


~ 


At high temperatures this reduces to 
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which, apartGirom the factor 2k,T, is identical with 
expression (4.1), recognizing of course, that the 
phonon frequencies are slightly different in the two 
cases. 

Ignoring this difference for present purposes, 


we see that Fay (a) is proportional to the angular 


averaged structure factor Say (a): 
Che) 
q 
Favea: < “ay ‘al = le Baa) tg: 


Greene and Kohn [32] have calculated this func- 
tion for Na and found that it was small and essentially 
constant for low q-values (normal processes) and rose 
steeply as one approached q = 2k, (umklapp processes). 
This behaviour is not peculiar to Na of course, and is 
a Characteristic feature of the structure factors. In 
general, one will find sharp peaks in S (q) whenever q 
is close to a reciprocal lattice vector, the shortest 
of which lies just above 2k, for the monovalent metals, 
Ane just below 2k, for, the polyvalent metals. 

Another feature of the results of [32] is that 
at high and low temperatures alike (273°K and 40°K), 
the ratio of the structure factor at gq = 2kp to Bee 


at q = 0 is unchanged, and in both cases is found to be: 
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We may readily obtain a qualitative understand- 


ing of these effects by considering the case when q 


~ 


is parallel to one of the reciprocal lattice vectors 


Gon irawe ignore dispersion, the vakue of (4.1) in this 


~ 


direction is given by 


Pe 
ied) > [oG- a) + ola - 9 i, 
2Msq (q-G) 


where Sc is the velocity of a longitudinal wave in the 
direction of G, and §@ is the unit step function. The 
£anction.in_curly,brackets.is;,equal to.unity up tothe 
zone boundary (q = G/2), and thereafter rises steeply 
aS gq approaches G. For Na (B.C.C. structure) the 
shortest G is the (110) direction, and Pete LLG of 


(4.2) evaluated at gq = 2kp to the value at q = 0 is 


given by 


with G110/7p = 1,148, “This result 1S the: same at’ all 
temperatures, since kp and G are both proportional to 
the inverse of the lattice constant. 

In order to do a sample calculation of? the 


electron-phonon (e.p.) enhancement factor for Na, we 
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Willeuse the functional) formigiven iby (4.2) in lace 
erytneraveraged)structure fLactor;smodified to the 


extent that we will define an effective G = G, such 


that 


-) = 40 


which gives G/2k,, e158, The behaviour of0(4. 2) 
forjthis value of 1G is; ildustrated in Api co enol CE BO 
atfbfters from the corrects form of the structure factor 


in that the latter is a little higher in the inter- 


mediate (q vu ky) region, this difference being mainly 


due to our neglect of dispersion. For polyvalent metals, 


OL #Gourse, Our approximation for Ss may not be used, 


AV 
because of the divergence at q = G, and one must doa 
more careful calculation in these cases. ‘Seé, for 
example, the resistivity calculations of Dynes and 


Carbotte [33]: 


Using (4.2) game (6. 5) and (3.16) (thus “gives 


Pee oe Ee 3 af2 
3) @H = SBE (Ey [fax x|w G0) |? 
Yo Gopi. F Ms 0 
I : : 
+ | ax Se | W(x) | | 
ae (x-a) 


where a = G/2kp, x = q/2k,, and n, (ER) = 3N/2E,- 
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Pig. léS Approximate structure’factor for Na. 


See equation 4.2. 
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A value for E./Ms* is obtained by calculating this 
quantity for longitudinal waves in the (100), (110) 
and (111) directions, and averaging the results. 
This gives a value for Na of E,,/Ms* = 1.22. 

For a pseudopotential, we will use the simple 


functional form suggested by Veljkovié and Slavié [34], 


namely 


(234), Banya LE 


which they fitted to the Heine-Abarenkov form factors 
[35,36] for a number of simple metals. In real space, 
this corresponds to a pseudopotential of the form: 


vir) = - (EB) WAS 6 (rer) 


defining a "pseudo-atom" of radius ee = b/2k,, which 
turns out to be a little smaller than the atomic radius. 
For Na, these authors give a value b = 27(.56121), and 


we obtain, upon performing the integrations in (4.3): 


(2h) = .20 (Na) 


with 303 of this value coming from the umklapp processes. 
This result compares well with the value .18 


given by Ashcroft and Wilkins [13], and if we add the 
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electron-electron contribution of .06 quoted by these 


authors, we arrive at 


(ley es, 9G (Na) 
Yo 


compared with the experimental value 1.25 quoted in 


[13]. 


§4.2 Binary Alloys - Interference and Impurity 


ConerrbDucions 


TO return towne binary substitutional ,alloy, 
we write the scattering function, with local pseudo- 


potentials, and treating the § C Dart wii ee nes eotatic. 


Cc 


approximation, as: 


Tiq,w) = |Wq)|* Soo (q,w) +c (1-e) [AW(q) [76 (w) 
where AW(q) = W, (q) - W, (q). We may note that since 
W = (l-c)W,+ cW,, we have 

2 oe 242 . 

W' = (l-c) W, +c Ww, + DoW c)W,W. 


ee (TSG (AWG 


2 2 
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and thus [(q,w) may be written as 
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When we use this expression to calculate corrections 
to the electronic.specific heat coefficient, we will 
have fourajeontributions.. The “frst two terms will 
give the linear interpolation between the electron- 
ion enhancement factors for the pure substances, while 
the third and fourth terms represent quadratic devia- 
tions due to dynamic and static effects respectively. 
It is interesting to note that the last two terms in 
(4.5) have the same concentration dependence, and 
involve the same combination of pseudopotentials. A 
Similar result is also true for v-component alloys. 
The contributions to the specific heat coeffi- 
cient due to the non-linear terms may be written, using 


PST Ley ana’ lost oO). aS 


2k 
| nee} F dQ Suny (G70) 
(AY) = -c(l-c) -2 | dq g|aw(a) [7/2] au a 
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aE 3 | aw(q) |* 
Ay. = Oo Feary “ad = 
Som. eh, 2 anv. dk [SUrameast: 
Yo imp Nk q70\Skedro. kK! kek, 
where int = "interference", imp = "impurity", and we 


have taken E, ok 7am: 

In order to investigate the qualitative behaviour 
and relative importance of these two terms, we will 
again use (4.2) as an approximation to the structure 


factor, and thus obtain, after some direct algebra: 


swoiveriod 6 eyidesttiap arts ptapltasivnt fo a obae! at 


Kae: aie | tasiisbat20s dno 4 oPE 


Lb emrse’ ond sortkt eat is 


-sivab i EGE amore: cnaes 94 
yiewhsosyeor etostisa sitete Dns: 


nk aneread ows teal Ssnt. jen edon, 


Bis ,Somebnedeb Ac vise stH88nS, ene 
A. ,elsitdetoyobuse, 20 coltatdamog ‘bated ‘ertonnt 
,eyor in sian ai lial tod ent oes’ ai’ “o fiheese ane ; 
~tt#60n teoc oft29 isqe edd oF -aotsudtediop- eam | 
pit ber, tee yng ed yan amyes te9nt f-nen ond lial i 


2m rad one) 


aw bas ("yi beieme” = qui \ "Batigas bait dai exeity 
nga = ‘ae nodes ved 


Ew: ‘BW! ante: awh oaatt ps ap fred 


ssvitbirsde etd 2) Aoktamiiouags 6: an. 6.8) oan akepe 
veadovts sods 5moe vette \atesale ed baie . Lede | i 


70 


an UE) 
ESO a ET (a Ghee) ene Be Ig 
Yo: int NEP ant 
2n 2) 
(AY) = [c(l-c) SHOOTER A - 
Yo imp NEP +P 
where we have defined: 
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It is interesting to note that the factor involv- 
ing the logarithms is the well known screening function 
for the electron gas. E,/Ms* in this context is to be 
viewed as an averaged of the pure metal values. 

The interference term is always negative, and thus 
constitutes a reduction of the electron-phonon enhance- 
ment upon alloying. The impurity term may be of either 
Sign ,) and tor |) AW (x)= constant, i is zero. “It will 
be positive if the integrand (apart from the factor 
ty (tT) is greater to the right of the singularity, 
anadunegative if it isngreater.toathe Pelt of ic. 

For illustrative purposes, we will examine a 
series of "theoretical" alloys in which each constituent 


has a pseudopotential of the form (4.4), thus giving for 
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sin 2785X sin 278 5x 
AW(x) = PL Saat pee Cr Pn eroneaxe: 4 


For most metals, the published information indi- 
cates that the first zero (q,) of the pseudopotentials 
occurs between gq = kn and gq = 2kp [7]. For the form 


which we have chosen, the first zero occurs at 


Thus, if we allow By and Bo to vary between 0.5 and 1.0, 
we will cover the range of interest. 

We will use the F.C.C. lattice structure, and 
in each case will calculate Ent for two values of "a", 
namely a = G344/2kp and a = Go 909/2kp: These are the 
two shortest vectors for this lattice, and should repre- 


sent upper and lower bounds for tent? 
Case l : By = Bo 


As a first example, we shall consider the case 


in which Qo is the same for both constituents, so that 


sin 2n8x)* 


2 
| ( 278X 


| AW (x) te = |AW(0) 
In this case, both nO and f imp are negative 
for all values of 8, and the results of the computa- 
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tions are given in Fig. 2, with E,/Ms ee he 
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The "correct" value for the interference term 
should lie somewhere between pees and cent subject, 
of course, to scaling, depending on whether E,,/Ms* is 
less «than, or greater than unity. 

in any case, we see that both the impurity and 
interference terms are of comparable magnitude, a con- 


clusion which will not likely be altered, in substance, 


for real materials. 


Case 2 ; W, (0) om W, (0) 


In this example, the pseudopotentials for each 
constituent have the same value at q = 0, but have 
different values of qo: The impurity contribution takes 
on both negative and positive values for different com- 
binations of By and Bo, asuinoicaled «in Big .43.. 

We see that for a fixed value of B, say, timp 
tends to negative values if By < B5 andowill thus act 
in the same direction as tent? For By > Bor however, 
the tendency is to positive values for fimp* One might 
thus expect some cancellation between f,,, and timp TOL 
these cases, since both functions are again in the same 
range as far as magnitude is concerned, as one may see 
from Fig. 4 and Fig. 5. These results are for 8, fixed 


at .65 and .85 respectively, with By ranging from) 0,5 
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One may also see that fin and Fimp tend to 


& 
assume their largest (absolute) values when By and 
B. are appreciably different from one another. This 
would be expected intuitively, since different values 
of these parameters are an indication of different 
scattering properties for the alloy constituents. 

The same trend is followed when W, (0) and W, (0) 


are allowed to vary at the same time as By and Bo, and 


the qualitative nature of the results are much the same. 


64.3. (AG-Au System — General censiderations 


Continuing from (4.5), we may write the elec- 


tronic specific heat coefficient as: 


ie ee alae) A, 


ine camp: 


G48.) yoo ace Co) hay 
where we have changed notation slightly, in order to 
write the concentration factors explicitly. 

If we now assume that the zero-order coefficient 
‘varies no more than linearly between the pure metals, 


and iwrite iLt.as: 
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Onsthe other hand, if we? collect all. the linear 


terms in (4.8), we may define: 


(4.10) ves (Lec) (yo) tay 2) + Pek altar 


ep fe) ep 
and it is seen that 
Pir wad 2) mig apie aes 
, ¥ (1=c) (yo say 2?) cy tay YY 
fe) fe) ep O ep 
Pore = 0 or | vs is equal to the measured value of 


the specific heat coefficients for the pure metals 
(including band structure and electron-phonon devia- 
tions from the free electron values), and is the linear 
interpolation between these values for c # 0 or l. 
Davis and Rayne [31] fitted their electronic 
specific heat measurements for Ag-Au alloys to an 


expression of the form 


(4912) ty Hey = K fob: c) 


in which ‘Veer was the linear interpolation between 
their experimental results for pure Ag and Au. A least 
squares fit to their data gave a value of K =-.138, 


with the largest deviations being at the Au end of the 
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alloy series, as one may see from Fig. 6, where we have 
plotted their results. 

These authors assumed Yo be) to be the zero-order 
“band structure" coefficient given by (4.9). This was 
consistent with the model of Stern [16], which they 
wished to test for concentration dependence, and in 
which the dynamic contributions are assumed to be zero. 

In point of fact, however, since Vdc) was 
obtained from the experimental data, it is clear that 
this function is identical with 15 as defined by (4.10), 
and already includes the linear electron-phonon enhance- 
ment. 

Thus, in reality, the function which Davis and 
Rayne fitted to their data is given by (4.11), and the 
coefficient "kK" determined by them contains the dynamic 
effects of the interaction through AY ser as well as 
the static contribution investigated by Stern,” and 
given in our case by AX ianp* We ee consider these 
effects more explicitly in the following sections. 

Note "that for cases in which the zero-order co- 
efficient has more than linear variation with concen- 
tration, this iwill- simply add another jtermto: ‘the 


numerator of (4.11), and may be thus readily included 


in the analysis. 
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-.138 c({l-c) 


Deviation from linearity (in concentration) 
of specific Neat coefficient for Ag-=-Au 
Points - Davis and Rayne experimental data. 
Solid: curve - least squares £1t to data. 


See equation 4.12. 
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84.4 Electron-Phonon Enhancement for Ag_ and Au 

FOrraAg and .Au, {!,C.iC.. structure), che first 
reciprocal lattice vector (Gi 44) is shorter than the 
shortest vector (Gi 40) for Na, and one could thus 
expect the structure factor to be relatively larger 
in the back-scattering region (q % 2k). We will use 
the form (4.2) to represent the structure factor, and 
in the absence of any published information, will 
choose a value of "a" such that Say (2Kp) /Say7 (0) is 
4/5 of the value calculated with the shortest reciprocal 
lattice vector, as in Na, This gives a = G/2k,= E123, 


and 
2k,- G 


Then, using expression (4.3), and the pseudopotentials 
Calculated: by Moriarty \137]. (Pic. +7), we_obtain, using 


E,,/Ms* =0.34 and 0.24 for Ag and Au respectively: 


(OX) one (Ag) 


Ove. Dp. 


(OY) vee (Au) 


103 GgDa 


The umklapp processes account for 77% of the value for 


Ag, and 86% for Au. 
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We have set ny (Ep) = (3N/2E,,) (m)/m), where my, 
is the effective mass that would be obtained from a 

band structure calculation. The existing band struc- 
ture calculations for Ag and Au are not in much 
agreement, and the published values for the band mass 

TM, (in units of m) vary between .88 and .98 for Ag; 

and, peeween .os and 1.01 for Au.” See [38,39,40] and 
references therein. 

If we use our calculated values for CON) ars and 
work backwards, we find band masses of .79 and .81 for 
Ag and Au respectively. These values seem somewhat 
low, of course, but in view of the existing disparity 
in the literature, one should use Moriarty's pseudo- 
potentials, temcalculate the. band structure as well, if 
one is to make a sensible comparison. As one may see 
from Fig. 7, Moriarty's pseudopotentials Sas quite large, 
especially in the backward direction, and may well alter 


the band mass (relative to the free electron value) more 


than the potentials used in previous calculations. 


§4.5, Interference ,Term for Ag-Au 
Following the same procedure as before, we use 
(4.7), with AW obtained from Moriarty's pseudopotentials, 


to calculate the interference contribution for the 
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Moriarty's pseudopotentials for the Ag-Au 
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Ag-Au system. We take E,/Ms* =mi29, ise sothe.simple 
average of the pure metal values, and thus arrive at 


=D 


s = -.037 c(1-c) (2) 


Yo int 
with this value coming entirely from the umklapp pro- 
cesses. This result is considerably smaller than it 
should be, and the problem appears to be with the 
pseudopotential differences used in the calculation. 
Indeed, if we use the same values for AW to calculate 


the residual resistivity of this alloy, we obtain: 


Sate 
BAe 97 hi c(1-c) | as x? | AW(x) |? 


2 2 
e (2mE.,) Ee 0 


= 9.65 c(l-c) pQ-cm 
in contrast with the experimental value 


a 35.06 c(l-c) pQ—-cm 


obtained by Davis and Rayne. Thus it would seem that 
the values of AW are too small, in the umklapp region 
at least, by almost a factor of two. This does not 
mean that the pseudopotentials for the pure metals are 
off by this much, of course, but since both of these 
functions rise steeply in the umklapp region, any small 


errors in each one will be magnified in the difference. 
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In fact, if one were to simply move the Ag pseudo- 


porentialvtouthe right "by 0,1 k the differences 


F! 
would increase by more than a factor of two in this 
region. 

Thus if we increase the values of AW by 
(35.06/9.65).2 ingthe calculation Of AY; nee we arrive 
at 


Mb 
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which is a substantial contribution to the experimen- 
tal value (-..138 c(1l-c)), even if the band mass is as 


lownasen 8. 


$4.6 Impurity Term for Ag-Au 


We Banner use AW from Moriarty's pseudopotentials 
to calculate AY imp! Since the values are not given for 
qs 2kp- We may note, however, from Fig. 7, that AW is 
practical ly zero for oa<, 1.5 kp and begins to increase, 
presumably to a peak just to the right of 2k At large 
values of gq the pseudopotentials will be again essen- 
tially zero. . Thus weymay tentatively expect: the “alloy- 


ing potential" (AW) for Ag-Au to be peaked at one value 


of g and to decrease rather quickly away from this value. 


There is some "circumstantial" evidence to sup- 


port this view if we consider the Cd-Mg-Hg system. All 
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three of these elements have essentially the same 
valence, size, and Fermi energy, just as do Ag and 
Au. “The data published by Inglesfield [41] on the 
alloys of these three metals demonstrates that the 
"alloying-wavenumber characteristic", which is 
essentially | aw| 2 times the screening function 
appearing in (4.7), has a pronounced peak for a 
particular value of q, and falls off rapidly other- 
wise. 

At any rate, to obtain some idea of what might 
be the case for Ag-Au, we take a simple functional 
form for | aw]? and calculate f,., and f,. We choose: 


ba imp 


2 
=(x/% 4} 
| AW(x) | 7= a’x* e ; 


which has a-single peak at gq = q1= 2k eX] 0 and perform 
the computations for qy varying between 0 and 4 Ker 
taking E,,/Ms* = .29, The principal value integration 


£Oxr en is. readily performed by noting that 


(or) Ph Ds fore 

== x a 2 x2 
| dx © 5 min ie Lm | e dx 
0 oe 0 


which is easily evaluated numerically. The results 

are. given in Fig: 8 as a function of qi: We see the 
characteristic: change otfisigny for fimp as the weight 

of the integrand shifts to the right, while an becomes 


monotonically more negative. 
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Fig. 8. Interference and impurity terms for the form 


of alloying potential considered in 84.6. 
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A reasonable inference from this calculation 
is that the impurity term for Ag-Au alloys is small, 
and possibly of either sign, and that the interference 
term causes most of the deviation from linearity 
(in concentration) of the specific heat coefficient 
ie Te=F0. 

At high temperatures, of course, only the static 
impurity term will survive, since the dynamic effects 
of the electron-ion interaction become negligible in 
this limit, as we have remarked earlier. This fact 
may, indeed, provide a means of experimentally deter- 
mining the relative importance of these two contribu- 
tions. The electron paramagnetic susceptibility, for 
instance, depends on the density of states, and can be 
measured both at high and low temperatures. At the 
present time, however, detailed data that would permit 


such a comparison do not seem to be available. 


SA.27 Concluding Remarks on Interference and Impurity 


Contributions 


From the results of the preceding sections, we 
may conclude that, in addition to the well-known impor- 
tance of the electron-phonon enhancement, the dynamic 
"interference" effects are as important as, and actually 


tend to dominate the static "impurity" effects at low 
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temperatures. This is mainly due to the fact that 
the umklapp processes give rise to large peaks in 
the structure factor close to the reciprocal lattice 
VECtCOrs. 
For monovalent metals, the shortest such vector 


is larger than.2k and thus only the steep "shoulder" 


FE? 
of,the.structure factor enters. theyscalculation«,,For 
the polyvalent metals, one or more of these peaks are 
included in the range of integration, and one may 
expect the dynamic effects to become even more pro- 
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CHAPTER 5 


SUMMARY AND CONCLUSIONS 


We have developed a formalism which describes, 
in second order perturbation theory, those thermody- 
namic properties of a metallic system arising from 
Ehe.electron-ion., interaction...The,centrak result is 
given by an expression for the free energy in terms 
of the electron-ion pseudopotentials and the dynamical 
structure factors of the ionic system. 

This result may be applied to the wide range of 
metals (and their solid and liquid alloys) which may 
be adequately treated by pseudopotential methods. By 
using appropriate dynamical structure factors (whether 
theoretically or experimentally determined) the forma- 
lism allows for arbitrary dynamic properties of the 
ionic system, and is valid for alloys in which mass 
differences and lattice distortion (due to volume 
differences etc.) play an important role. 

A significant aspect of the central formula is 
that it embodies, in a single concise expression, 
many results which are usually treated quite separately, 
using a variety of different techniques, which are 


often difficult to extend without a complete reformu- 


lation. 
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We have shown how one may simply obtain a 
number of such results from the general formalism, 
and have further derived a concise expression for 
the low-temperature electronic specific heat changes 
in a substitutional alloy, due to the dynamic part 
of the electron-ion interaction. In the numerical 
calculations we use this formula to perform some 
Simple illustrative calculations of the electron- 
phonon enhancement factor. 

An interesting part of the calculations was 
concerned with a comparison of the relative impor- 
tance of the dynamic "interference" and static 
"impurity" contributions to the low-temperature 
electronic specific heat coefficient in binary alloys. 

For a series of "theoretical" alloys, both of 
these terms were found to fall generally within the 
same order of magnitude, with the interference term 
demonstrating a tendency to dominate. The importance 
of this term is largely due to the. fact that the 
nen processes give rise to pronounced peaks in 
the structure factors near reciprocal lattice vectors, 
an effect which is well known in resitivity calcula- 
tions. 

For Ag-Au alloys, the deviation from linearity 


(in concentration) of the low-temperature electronic 
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specific heat coefficient was examined, and our results 
seem to indicate that the interference term in this 
case is capable of adequately accounting for both the 
general magnitude and concentration dependence of the 
observed behaviour. 

The relative importance of the interference and 
impurity terms might actually be determined experimen- 
tally from electron paramagnetic susceptibility mea- 
surements. This property is proportional to the 
density of states, and is measurable both at high 
and low temperatures. At high temperatures the inter- 
ference term, along with the other dynamic parts of 
the electron-ion interaction, become negligible, while 
the static impurity term remains essentially unchanged. 
One might thus be able to separate the two terms in 


this way, if detailed experimental data were available. 


Se 


atitest wo Baw benimass any snatoiMaeeg 2 ' 
aid? ot med Sone t9Ixedake of tery ON 2 

oft tod 202 pabinvooos eloseepebe’ 20. @idagno et | a i 
att 26 epasbhbnegsb notiswolssaew inp abit ingsm serenee 
bas sunszoiteini ai3 iq somatiogns ovianten sit a 
~nesrizegxne hentmrez'eb eel PRlevIos tipim aarams er 
~san yiilidisqercese »isenpemeredg nostne te mers t ‘yited 
ed 03 Letoiduoscng é) yaeeQete aaae -edaomoine | 

| ipid te dtod oldsxvesom &f bas “(aeaase 30 oan 
~xederd ) Stic esruterounet dein 2A haat niet el wok bite 
%© ettag a ain ah reiito ens ntatiw paalea ilies eonese2 " . 
aL.trlw \ekdivélpen anoved *, tokiugionak noi~nentaete wet | a 
-bopusonm ylisttneare ankamey mee ytkawgmk obdadey “7 i" 
mi eames owt sdt sieuayee of Side! od euds Jitptm en 
,oldslieve siew sapeb cotenalnaiaie befistab ti wesw eid 


mike ¢ 


ti) 


[2] 
gel 
[4] 
[5] 


[6] 
[7] 


[8] 
[9] 
[10] 


[11] 
ee 


13) 


[14] 


REFERENCES 


N.F. Mott and H. Jones, "The Theory of the Proper- 
ties of Metals and Alloys", Oxford University 
Press, London (1936). 

Ee BLOCH, 2.) nVYork 92, 999 1928) « 

Cy. Herring, Euys. KEV. 27, LLOd. (1940). 

Cy Herrrng anc A.G.e Hill, Phys. Rev. oo, 132 (1940). 

J.C. Phillips and L. Kleinman, Phys. Rev. -116, 287 
CES So). 

He anconcik , J. Pnys. Chem.  oOlLLdS LO; RSI aeKLI5 9): 

M.L. Cohen, V. Heine, D. Weaire, "Adv. in Solid 
State Physics’, ed.) HH. Ehrenreich, F. Seitz, 

D. Turnbull, Academic Press, New York, Vol.24 
RIO es 

Re Penerls, sae yolk C0ye 1 Osen(Lo35):. 

L. Van Hove, Phys. Rev. 95, 249 (1954). 

Mio, puckwingnam andr M. Re ochalroth, «Proc. PNYS. sOC. 
AO ao 2o Wy oo 4)r. 

RB. Migdalyy sovaetePhyserJELP Jy) 990" (1958)". 

S. Nakajima and M. Watabe, Prog. Theor. Phys. 30, 
TIL 19 6S) eer LOGge i Neou. PHYS wat, woe. (1963). 

NeW. AsnerortrandJ.W. Wilkins, ePhys.» Lett... 147 
Zoo VL GS: ) is 


Povey rice Alms Physi 31, 100) (1965). 


93 


~reqoxe sta, 2q oni ort 28not : candied 
yttenewv ia bretxO ‘ “ aehLA, Bint ale av a 

ie teen), ob, 

»(8seL) eee ise “thewda: yt 

oneny: Call \d Jwen se eoizeet at 

COME) SEL, B2 von seyda \ DLs DLeien pe ij & 
8S abl wwSa Jayna ene BA, =k bata y euacidy ane a 

| ie wlOee sh oun a 

(CAL) BLE OL aberge me YAS .%; © aionosmn it! MB1 | 

hit tod i ¥bA" estes of SESH eV yoedoD yaa con si 
,stios ‘% lo herbseele «fh «Be Seakeyad aoegee eye | 

bs, Lav or wall ,aeer7". ‘eMehAoA, sh eaterea ae a - 

. (ERE) cay oe ales & valaohet th 81 

(ReOLy Spe ae VOR . ext oo as ak LET 

008 .aydt .oond iat uM aia medpagsoug +.M - 
| on (BORE) BRO haa: | om 
(8808) age 8 ateD ‘aydd dove’ , IebpdM .aa ign os 

GE .@yHt . 98d .pord ae 3H ‘M bus smitetel .2 ist) vg > 
(SRK), Le ‘= Bit: haan’ poms tenes) sty _ 


_ oF 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[27] 


[28] 


[29] 


94 


J.W. Wilkins, "Observable Many-Body Effects in 
Metals", Nordita, Copenhagen (1968). 

E.A. Stern, Phys. Rev. 144, 545 (1966). 

E. Haga yiProes Phys. rSoce ,9DjiyhG9 461967). 

M.R.«Schafroth, Helv. Phys. Acta. 24, 645 (1951). 

D.J. Thouless, "The Quantum Mechanics of Many-Body 
Systems", Academic Press, New York (1972). 

A.B. Bhatia and W.P. O'Leary, Lett.al Nuovo Cimento 
3 pelagic 72 a: 

A.B. Bhatia and D.E. Thornton, Phys. Rev. B2, 3004 
(ESROde. 

A.B. Bhatia (private communication). 

A, Bijieiatia and D.E. Thornton, Phys. Rev. B4, 2325 
619 Tia 

Ke, Huang, Proce Rey, £ Soet (At90 724102) . GL9 457). 

M.AsKrivogilazy "Theory iofiX\-Ray and Treen Neutron 
Scattering by Real Crystals", Plenum, New York 
CL9G Se. 

W. Hume-Rothery, G.W. Mabbott, K.M. Channel-Evans, 
Phit.) Thavs.Roy. Soc. .255A7, 1) (1954). 

Ce. ‘Kittel. "Ouantum Theory cf Solids, iWwiley, 
New York (1963). 

A. Messiah, "Quantum Mechanics", Wiley, New York 
(196 6) % 

be Sham and 0 .Mwaszaman, “Adv. in Solna State 
Physics", .ed.. /. seitz and) D. Turnbull, Academic 


Press, Newrvork, Vol. 25: (1963). 


"i" ye ag a b Ja = PDT ne ‘ Sale ye 
iF +) fi 7 7 , ' 1 . aw _ 
7 Wy AY yi hoe | ane “_e : ; 
; . 6 OS @ Toe _ 


is aaa, 1% iy 


nt aseehas sioe~ erat ofisvasede" | 
(8820) popbiueies ‘wl 


ag 1 inet, a ae 
(1802) 20 ohS pecan -Byds wed Asedandoe AM 
yboa-ynsM ‘to eo inet: 79M mudnsyd enh” .aekbodit bed ro) | 
. (ST OL) #20Y walt \2@e it LMP GON . tad a | me 


.(Wa@L) GAL (Fe .068) Lauda Lois 


otaemi oven Is.j3tad.,yxsed'O «ult bus pbctienlls os tos 
| CPO Me 9b | 

SOE , SH. .von-.2ute othaod® ted bey Bitade, GA wee mi 

PAMPER ae 

(tolysotnummios sdevtxe) shin Buh iss} 

CSEg , 8a ven dyna ionized? (Gs. bene. siseda. AK (esi 

: (AV @L) ‘ Wy 7 

~(ThOL) SOF ORta* ode Yee 208 ema zt ts) 

coxjue4 LemzeiT bas Yaak to yeoman"  SOLROW EAH AM (2s) c | 


ae ) Hi ie ay 
oh 


| Seeaenp oy rt 

cadevE~Lonisio «MM. it idaddtast We ] : adaelily we —— 4 
-CPEGE) L ARES 302! You eng Liat “ah a : 

syed iW Nabi to® t6 {root meu” iatt Et ” ray Oh ‘ 

| ays (BBR) Haot wom ete 
shenviad weit a "sotnedoen it ei . 4 $} ia ; 


wy or - 


- 


AxoOY wot ny mobs ,"alssieyt> Leek ed el 


’ 


es 
Sgn7@ Bited nt <¥BA" vnsmiEs) MG ctw = ea 
anaes ce cal Hite ayia oY a + teat 


) — 


[30] 


[31] 


[32] 


[33] 


[34] 


P25) 


[36] 


Rea 


[38] 


[39] 
[40] 


[42] 


20 


NiF.- Mott and E.A. Davis, “Electronic Processés 
in Non-Crystalline Materials", Oxford University 
Press, London (1971). 

THe Days |ang.J.A. Rayne, Phys. Rev. B6, 2931 
72) vs. 

MoD Greene and W; Kohn, Phys. Rev. 137), 513° (1965). 

RC. (bynes ana J.P). Carbotte, Phys. Rev. 175; 913 
(1969). 

We Weljroviceand Slavic, Phys. Rev. .Lett. 29, 105 
(E972). 

Wo Heine and Eb. “Abarenkov, Phil. Mag. 9, 451, (1964). 

A.O.E. Animalu and V. Heine, Phil. Mag. 12, 1249 
CLO 6: 5 )F 

J.As wondarty ehyee Rey, (Bl) 1362. (1969). 

Wed. O' SulMiven, A.C. Sswitendick, JsE. Schirber, 
Phys .. Rev.) Bl, 144301970)". | 

PE. Lewiswand, P.M. bee, Phys. Rev. 1/5; 795. (1968). 

NBs  Cherstencenymphyc. Stat. Sol. (hb) 54," 551 (1972). 


J.B. Ingleéstield, Acta. Metall. 17,1395 (1969);. 


_! 


ne 


Men iy 


me ay 


a 
Ps 0 - ny 1g 
ry yy im ie aN i , 
. - ei he aan age 
1 
1 
i 
1 
bil 
ae 
i 
i L 7 7 ry / 
nn 7 
; ’ 
' 
AY, i] : mS i 
Pe 
, > 
y i 
TIA aoe n 
f \ 


at ‘ai! on 
tp mt , 


